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Abstract. Numerical studies on the unjamming packing fraction of bi- and polydisperse disk packings, 
which are generated through compression of a monodisperse crystal, are presented. In bidisperse systems, 
a fraction /+ = 0.400 up to 0.800 of the total number of particles have their radii increased by AR, while 
the rest has their radii decreased by the same amount. Polydisperse packings are prepared by changing all 
particle radii according to a uniform distribution in the range [— AR, AR]. The results indicate that the 
critical packing fraction is never larger than the value for the initial monodisperse crystal, <j>o — 7r/vl2, 
and that the lowest value achieved is approximately the one for random close packing. These results are 
seen as a consequence of the interplay between the increase in small-small particle contacts and the local 
crystalline order provided by the large-large particle contacts. 

PACS. 05.10.-a Computational methods in statistical physics and nonlinear dynamics - 64.70.kj Glasses 
- 64. 70.ps Granules 

1 Introduction or coagulated colloidal microstructures with high elastic 

moduli are a few examples with practical applications that 
The jammed state of condensed matter is characterized by exllibit j ammin S- ° n the theoretical side, random pack- 
the sudden arrest of a system's internal dynamics. Macro- in § s of hard and soft elem ents (spheres, disks, etc.) display 
scopically, the jammed system develops an yield stress and similar behavior to those observed in real systems and are 
behaves, essentially, like a solid. A pile of sand under the often used as Prototype systems for the study of jamming, 
action of gravity, clogged flow of powders through a pipe lt has been extensively studied recently both by Simula- 
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tions [UHHIlElEJIZllHlEEn] and experiments QTJQ1]. As 
expected for strongly interacting systems, there are sev- 
eral theoretical proposals [T5l[T^l [T5 , 16 to describe this 
state, most of them relying on approximate mean-field ar- 
guments. 

There are some facts that seem well settled today about 
jamming. First, a quench is an essential ingredient for a 
system to reach the jammed state. It prevents any crys- 
tallization that may occur during a slow rearrangement of 
particle positions during equilibration. Second, the criti- 
cal jamming density is affected by particle size ratio |10j . 
shape [17] and the preparation protocol, but its critical 
properties are the same (recently, Chaudhuri et al. [9] 
showed that this critical point is not unique, even for large 
systems). Third, the jamming point in monodisperse and 
bidisperse packings is manifested in structural properties 
as the 5-function behavior of the first peak of the radial 
distribution function and the split of its second neighbor 
peak [7J. 

Most studies about jamming at zero temperature fo- 
cus only in one preparation protocol (random packings of, 
possibly, overlapping particles, quenched to the nearest 
minimum energy state). One may ask, then, is it possi- 
ble to produce a jammed state from a completely ordered 
system (crystal)? Is the quench alone enough to take the 
packing out of its global minimum of energy and trap it 
in a jammed state? In particular, if the jammed system 
prepared in this way will it be more dense than the ini- 
tial crystal? All these question, and others, will be ad- 
dressed in this paper. Jammed states will be produced by 
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the quenching, and further decompression, of crystalline 
disk packings. 

This initial condition, at first sight, is not suitable to 
produce a jammed state, since disordered packings are 
commonly associated with jamming. However, it will be 
shown that the structural features of jamming are present 
in such systems, hence regarding this preparation proto- 
col as a valid one to produce jammed packings. It will 
be shown that this initial condition opens the possibility 
to reach jammed states in a distinct region of the pack- 
ing phase diagram, which are unaccessible from ordinary, 
random initial packings algorithms. Along with the studies 
of bidisperse packings, the decompression of polydisperse 
packings is explored. 

Sect. [5] is where simulation details are provided, sect. 
[5] holds all results and sect. 2] is reserved for conclusion. 

2 Simulation methods 

The particles are soft, elastic disks, which interact through 
linear springs. The compression potential energy between 
disks i and j given by: 

U tj = ^K(Ri + R 3 - r ij f9{R l + R j - m), (1) 

where, k is the elastic constant (taken k — 1 and equal for 
all contacts), Ri is the i-th disk radius, r,-j is the distance 
between the disks, and &(x) is the step function. 

Initially, the packing consists in N disks, of radius 
Ro (taken as the unit length), arranged in a triangular 
lattice. The contact energy, eq. ([IJ, is given in units of 
kRq. The system's periodic boundary lengths are given by 
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L x = 2R N X and L Y = V3R N Y , where N x x N Y = N, 
is the total number of disks. The values N x = Ny = 50 
are chosen, which gives N = 2500, throughout the ex- 
periments. The results presented here are essentially the 
same for systems with N x = Ny = 10 and N = 100. 
This choice of boundary lengths perfectly accommodates 
a triangular lattice of equal disks, which implies that the 
initial packing fraction has the largest value for a two di- 
mensional monodisperse system (the Kepler conjecture) 
US]: 

The quench is performed by changing the disks radii 
by a suitable amount, the dispersity degree AR. Such 
change is instantaneous, in order to trap the system in 
a jammed state (quench). For bidisperse packings, a num- 
ber N + = f+N of the disks (randomly chosen) have their 
radii increased by AR while the rest of the particles, 7V_ = 
(1 — ,f+)N, have their radii decreased by AR (similar to 
what was used in \19\). The number fractions chosen are 
/+ = 0.400, 0.500, 0.600, 0.700 and 0.800. For polydisperse 
packings, the radii are changed by an amount uniformly 
distributed between [—AR, AR]. Fig. [T] shows a quenched 
configuration. 

These changes introduce a compression potential en- 
ergy, since there will always be some overlap between near- 
est neighbor grown disks. Given the absence of energy 
dissipation, the minimization is performed by the conju- 
gate gradient method [20] , A dissipative packing cannot 
be studied with this numerical method. A molecular dy- 
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Fig. 1. An illustration of a quenched packing (/+ = 0.500 and 
AR/Ro = 0.300). 

namics (MD) approach is more suitable, and certainly will 
yield very distinct results (for a study on the phase di- 
agram of dissipative packings see |21|). Finally, random 
packings can also be prepared with a MD approach by 
swelling void particles (void expansion method [32]) that 
increase the volume fraction and takes the system through 
the jammed point. 

It should noticed that a jammed state is not reached, 
with this protocol, at lower number fraction, for instance 
/+ < 0.300. Given the small probability to form large- 
large contacts at the beginning, the available space for re- 
arrangements is larger than the one needed for complete 
relaxation. This leads to a melting of the packing, instead 
of jamming (the highlight emphasizes that the melting 
picture should be supported by additional simulations). 
Therefore, such low number fraction packings do not be- 
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have as their high /+ counterparts. This will happen only 
if the mean radius of such packings is larger than Rq. 

Since the goal is to find the maximum packing fraction 
with a vanishing potential energy, an initial minimization 
is performed right after the quench. The energy minimum 
is achieved when the difference between the current and 
the last energy values is no more than 10 -10 . After reach- 
ing the nearest energy minimum, particles are slowly de- 
compressed, i.e., particle's radii are decreased by a small 
constant amount 7 at each cycle, which provides a slightly 
larger space for them to relax (expansion step) . After each 
expansion step, an energy minimization is performed in or- 
der to take the system closer to the zero potential energy 
state. The decompression is finished when the total poten- 
tial energy is less than a predefined value, e. Hence, at the 
of the protocol for a suitable e value, the system should be 
very close to the jamming point. In sect. 13. ll the influence 
of both parameters on the results will be shown. 

The average packing fraction (</>), from now on referred 
to as the critical packing fraction (CPF), is measured at 
the end of the decompression as a function of AR/Rq. The 
averages are over realizations (typically 20 for each case) 
and the error bars are calculated as \J {{<f) ~ {(p)) 2 ) in each 
case. Also, the packing structure and order are studied 
through the calculation of the Radial Distribution Func- 
tion (RDF), g(r), and the orientational order parameter: 

Zj 

^ = —J2e m ' k , (3) 
Zj fe=i 

where the sum runs over the Zj nearest neighbors of disk j 
and 6jk is the angle between the line joining the j-th and 
fc-th disks centers and the x axis [33] ■ The absolute value of 
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tyj is measured at the end of the full minimization, and is 
presented as an average over particles and runs. Its value 
is unity for a perfect triangular array of particle, while 
< 1 for disordered packings. Two disks are considered 
first neighbors if they overlap. 

It should be noticed that, in most studies of bidisperse 
packings 2 ,3, 5 ,9,10 , particle size differences are given in 
terms of the size ratio, a, instead of the size difference, 
AR. Both quantities are connected by: 

AR _ a - 1 
Ro ~ o + 1 ' 

3 Results and discussion 

This sect, holds all numerical results. First, some results 
for the CPF as a function of the simulation parameters 7 
and e are shown. Second, the full CPF results, along with 
the packing structure and order will follow, respectively. 

3.1 Independence on the simulation parameters 

The simulation is controlled by two sensitive parameters, 
the decompression rate, 7, and the minimum compression 
energy, e. Figure [D holds the results for the CPF for three 
distinct parameter sets: 

The three cases presented, all for / + = 0.500, have the 
following set of parameter values: case 1, 7 = 10~ 5 and 
e = 10~ 6 ; case 2, 7 = 10~ 6 and e = 10~ 6 ; case 3, 7 = lO" 6 
and e = 10~ 8 . As seen in this figure, all results agree well 
within simulation error. Hence, all the following results 
will be given for case 1, unless noticed otherwise. This 
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Fig. 2. CPF for distinct values of the decompression parame- Fig. 3. All results for the CPF, for distinct particle number 



ter, 7, and the energy minimum, e. 

value of e corresponds to a minimum energy per particle 
of the order 10~ 9 . Similar results hold for the polydisperse 
packing (not shown). 

3.2 Jamming packing fraction 

In fig. [3l the results for the CPF for all number fraction, 
/+, and size dispersity, ARj ' Rq, values are shown. 

From this graph, one can see that the final packing 
fraction is never larger than the one for the triangular 
lattice, eq. and it goes through a minimum. From 
the inset of fig. |3l it can be seen that the minimum CPF 
changes with number fraction. Its lowest value is 0.843 
for /+ = 0.400 at AR/R = 0.120, while the largest mini- 
mum is 0.845 for /+ = 0.800 at AR/R = 0.150. The CPF 
value for the random close packing (RCP) state is shown 
as a reference value, since it is the jamming packing frac- 
tion for the monodisperse packing [2J. These results indi- 
cate that the RCP state corresponds to the lowest critical 
packing fraction achieved with this protocol, with packing 



fraction, /+, and size dispersity. Symbols correspond to number 
fractions of /+ = 0.400 (circles), /+ = 0.500 (squares), /+ = 
0.600 (diamonds), /+ = 0.700 (triangles), /+ = 0.800 (left 
triangles) and polydisperse packing (inverted triangles). The 
inset is a zoom to the region close to the curves minima and 
the dashed line represents the CPF value for the RCP state 
0, 4>rcp = 0.842. 

properties distinct from the original study (in that case 
2 j, /+ = 0.500 and AR/Rq = 1/6). For the polydisperse 
packing, the minimum CPF is 0.844 at AR/R a = 0.200. 

Recent results on jamming of bidisperse sphere pack- 
ings [I]]] and ellipsoid and dimer packings [17] . that fo- 
cused on the influence of the size ratio in the CPF, show 
that this quantity presents a maximum in the size ratio 
range [1, oo]. These references also carried out their simu- 
lations with random initial packings. (In [10] . the largest 
CPF occurs at /+ = 0.500 and AR/R Q = 1/3). Here, 
all CPF values are also above the RCP one, but it goes 
through a minimum instead of a maximum. The reason 
behind this distinction seems to be mainly the initial pack- 
ing. The jamming point in a random monodisperse pack- 
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ing may be increased with an appropriate size ratio (the The value of this mean radius is R — Rq — r, with r — jn, 
system packs more efficiently). On the other hand, the where n is the number of decompression steps. Defining 
jamming point for a regular packing can only be decreased x — AR/Rq and y — r/Ro, and dividing both sides by 
from the value given in eq. ([2])- This can be seen as a con- LxLy one has: 

A. 

sequence of the fact that, as stated earlier, the triangular 
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lattice is the most dense packing of equal disks [To]- There- Uji q _ jn) 2 + 2(2/ — l)(i?o — jn)AR + AR 2 ] 

LxLy 



fore, one can conclude that performing a decompression 
simulation with a regular initial packing, one can reach a 
very dense jammed state, not accessible from a random 
initial packing. This is the main results of this paper and 
it is what is meant by a novel route to unjamming in the 
title. 

The small dispersity behavior of the CPF can be un- 
derstood as follows. When the regular triangular array of 
disks is quenched, the total area occupied by the disks, 
initially given by Aq — NttRq, is changed to: 

' N + N- 



A b0 = n 



A 



ovlp • 



where A ov i p represents the total overlapped area. Develop- 
ing the squared terms and using the fact that /+ + /_ = 1, 
one reaches the following expression for the initial modi- 
fied area: 



The final average CPF, among different experiments with 
fixed boundary lengths can be obtained using the average 
r value in Ab- Therefore, using eq. ©, this relationship 
yields: 

(&> = 

0o [1 + 2(2/+ - l)x + x 2 - 2(2/+ - l)x (y) -2(y) + (y 2 )] 

(4) 

A similar argument holds for a polydisperse packing. 
In this case, the area occupied by the disks after the initial 
compression is given by: 

JV 

ApQ = 7T J2( R + AR if - A Ovlp, 
8=1 

where ARi is the change in the i-th disk radius. At the 
end of the decompression phase, the total disk area is: 



A p = irY,( R + AR i? 



A b0 = NirRl + 2(2/+ - l)NirR AR + NirAR 2 - A ovlp . 



with R given as above. Since the quantity ARi is uni- 
formly distributed in the range [—AR, AR], one may use 
Since the experiment is carried through particle decom- the moments of ARi, 



pression at a constant rate, when the overlapped area van- i 

_$>^ = 0, (5) 

ishes, the total area occupied by the disks should be given 
by the exact same expression, but for a distinct mean ra- 
dius R, instead of Rq: 



N 



i=l 



to write 

A b = NnR 2 + 2(2/+ - l)NnRAR + NnAR 2 . A p = N * 



A f? 2 

(i? -7«) 2 + ^ 
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Following the same steps as in the bidisperse case, the 
average CPF can be written as: 



t>p) = <Pa 



l-2(y) + (y 2 ) + 



(7) 



Hence, with the knowledge of the average number of de- 
compression steps, one can match equations (j4]) and ([7|l 
with the results given in fig. [3] The results for (y) and 
(y 2 ) are given in fig. 01 
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Fig. 4. Upper panel: average number of decompression cy- 
cles. Symbols follows convention in fig. [3] Lower panel: average 
square number of decompression cycles. 

Consider the following argument to explain this re- 
sult. If this decompression experiment was performed in 
the monodisperse case, where each disk has exactly 6 con- 
tacts, simetrically placed around its center, no change in 
the structure will ever occur due to the decompression, 
and the packing fraction value in which the compression 
energy is zero would be given by @. Therefore, the disk 
radii should return to their original value in order to reach 
this packing fraction. This implies that (y 7 ) = x 1 . Since at 



Table 1. Decompression step parameters, measured from the 
curves in fig. [4] 



/+ 


0.400 


0.500 


0.600 


0.700 


0.800 


poly 


a 


0.214 


0.413 


0.588 


0.722 


0.895 


0.244 


b 


0.046 


0.171 


0.346 


0.520 


0.802 


0.061 


a 


0.898 


0.941 


0.956 


0.956 


0.979 


0.941 


P 


1.80 


1.88 


1.91 


1.91 


1.96 


1.89 



low dispersity this relationship is approximately realized, 
one can infer that at the jamming point, one has (y) = ax a 
and (y 2 ) = bx^ . These four parameters a, b, a, (3 represent 
the effect of structure rearrangements in the quantity r. 
Table [T] holds their values as measured from power law fits 
to the curves in fig. 0J Since deviation from power law be- 
havior do not occur at the same dispersity for all cases, the 
fits were performed up to x = 0.040 for /+ = 0.400, 0.050 
for /+ = 0.500, 0.060 for /+ = 0.600 and x = 0.100 for the 
other cases. The parameters approach their monodisperse 
values as one increases the number fraction. The polydis- 
perse packing is an exception since it has no monodisperse 
limit. Also, one can see that, at all number fractions and 
the polydisperse case, the relationship j3 = 2a holds. The 
data do not imply a simple relationship between the coef- 
ficients a and b. 

By using this form for the average number of decom- 
pression steps, the average CPF, eq. igj), can be cast in 
the following form: 



^ = l+2{2f + -l)x+x 2 -2{2f+-l)ax a+1 -2ax a +bx 2a . 



(8) 
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This curve for the parameters a, b, a and (5 given for the 
/_!_ = 0.700 case is shown as a dashed line in fig. [5] One 
can see that, at low dispersity, this eq. agrees well with 
the results. 




Fig. 5. Eq. (J8]) plotted with the parameters correspondent to 
the /+ = 0.700 case. 

For the polydisperse packing, the average CPF can be 
written as: 

= i-2ax a + bx 2a + -x 2 . (9) 
9o 3 

The comparison between figs. [3] and |4] shows that the 
deviation from the power law behavior of (y) coincides 
with the region where the CPF reaches its minimum value. 
Since a complete understanding of this fact would be given 
by a more detailed (and complicated) theoretical approach, 
such as the one proposed in [16] . only the packing struc- 
ture is probed here. The reason behind this choice is that 
the structure relaxation during decompression is the main 
cause for the results seen in figs. [3] and |U 
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3.3 Packing structure 

The RDF is measured at the end of the full minimiza- 
tion process. These measurements are performed regard- 
ing the type of particle contact, i.e., the probabilities of 
small-small, gss( r ), large-large, <?ll(V), and small-large, 
9sl{t), contacts are measured. In the triangular array of 
monodisperse disks, one expects g(r) to have sharp peaks 
at r = 1, V3 and 2 particle diameters. Figure |H] shows all 
RDFs for AR/Ro = 0.120 and all number fractions. Fig- 
ure[7]has these same functions but for a dispersity value of 
0.500. In both cases, all interparticle distances are normal- 
ized by the corresponding final average particle diameter. 
For instance, for a small-small contact, the final small par- 
ticle diameter is 2 (R) = 2(Rq — AR — (r)), where (r) is the 
average number of decompression steps, fig. 2] Hence, the 
distance range where these contact probabilities are mea- 
sured appear distinct for each contact type and number 
fraction, especially for large-large contacts, which have a 
larger mean particle diameter. 

In all curves in fig- El the jamming structural signature, 
a delta-like first peak and a split second peak, are seen. 
Also, the second small-small contact peaks are not located 
precisely at r = \/3 and 2. Instead, they are shifted to the 
right, consistent with the results shown in 10 . A brief 
explanation of this fact is that these two peaks occur only 
when three particles form a triangular cluster. Hence, this 
type of cluster should be absent for small-small and small- 
large contacts. Moreover, small-large contacts have a sec- 
ond peak at intermediate positions between small-small 
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Fig. 6. RDFs for a dispersity value of 0.120. Each line, shifted 
for clarity, corresponds to a distinct number fraction (from top 
to bottom): /+ = 0.800, 0.700, 0.600, 0.500, 0.400. 

and large-large contacts, as expected, since the mean di- 
ameter of such a contact is given by 2(Rq — (r)). 
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Fig. 7. RDFs for a dispersity value of 0.500. Curves are shown 
in accordance to the convention in fig. [5] 



A marked feature of these graphs is that the gss{i") 
and gsL(f) peaks between r = 1 and decrease with 
/ + , while those for r — 2 become sharper. This indicates 
that, for larger / + , small particle relax to positions farther 



away from each other, even though they start all at the 
same structure. 

The <7ii(r) peaks at r = 1, and 2, increase and be- 
come sharper, while the intermediate ones between r = 1 
and \/3 almost disappear at the largest number fraction. 
Since sharp peaks at y3 and 2 are a signature of a triangu- 
lar lattice structure, one may infer that such large particle 
rich packings relax to structures similar to a crystalline 
one. This fact also explains why the large-large peaks be- 
tween r = 1 and y/3 become smaller. Such characteristic 
should be expected, since the initial packing is regular 
and only large-large particle contacts, at the outset, im- 
ply compression. Therefore, if a group of first neighbors 
become large particles, they will probably remain in this 
cluster up to the end of the process. 

An illustration of such a configuration is given in fig. [8] 
It shows a well mixed packing with an occasional pocket 
of large particle crystals (lower left corner) . 

On the other hand, the pair correlation functions at 
large dispersity, fig. shows markedly distinct features 
of the packing structure. First of all, small particle ag- 
gregates become progressively more distant, in small par- 
ticle mean diameter units, for larger / + . Second, 
at /+ = 0.400 and / + = 0.500 shows several small peaks 
between those at r = 1 and 2, as in fig. O but the one 
at r — \/3 is not easy to distinguish. Only at high num- 
ber fractions this peak becomes clear. This means that 
large particles form, again, structures close to crystalline 
ones at high number fraction. Finally, one can see a clear 
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Fig. 8. Packing configuration for AR/Ro = 0.12 and /+ = Fig. 9. Packing configuration for AR/Ro = 0.500 and /+ = 
0.50. Particle periodic images are omitted. 0.800. Small particles are enlarged for better visualization. 



change in gsLij 1 ) from / + = 0.700 to 0.800. The peaks 
for the lower number fraction are sharper than the corre- 
sponding ones at larger number fraction. This is probably 
due to the fact that small particles can be more easily 
accommodated in vacancies among the large particle con- 
tact network (approximately triangular) . Figure [5] shows a 
snapshot of a packing at the largest number fraction and 
dispersity. One clearly sees that the structure is a crystal 
with defects. Finally, the broad gss( r ) peak around r = 1, 
at / + = 0.800, does not imply overlap between small par- 
ticles. In fact, the gss( r ) value is zero while r < 1. The 
reason for this apparent broad peak is that the bin size is 
larger for smaller average contact diameter. This quantity 
decreases for increasing number fraction and dispersity. 
Hence, the bin size at /+ = 0.800 is significantly larger 
than those at lower number fractions. 



Although not seen in the RDF plots, the area below the 
first peak for all RDFs changes as a function of the number 
fraction and dispersity. This might have implications for 
the mechanism that leads to the CPF values seen in fig. 
131 In order to study this influence, the area below the first 
peak, in each of the RDFs, is measured. The r range in 
which this area is considered is from the first peak position 
(contact diameter) up to this distance plus the dispersity, 
AR. This is an (unnormalized) account for the average 
number of neighbors of a given type around a given parti- 
cle [7] . A more natural approach to this question would be 
to compare the coordination numbers with regard to each 
contact type. Any protocol for producing jammed states is 
known to produce an amount of rattlers (particles with no 
contacts). These particles also contribute to the packing 
relaxation. Then, a comparison of coordination number 
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would exclude these particles from the analysis. The no- 
tation corresponds to 

dss+AR 

Nss(AR) = J 9ss(r)dr, 

dss 
d LL +AR 

N LL {AR) = J g LL (r)dr, 

dLL 

dsL+AR 

Nsl(AR) = J g S s(r)dr. 

dsL 

This choice for the integration limits ensures that, for 
small-small contacts, only small particles are within this 
range. At the first few dispersity values, which are of the 
order of the g(r) bin size, this integration leads to an over- 
estimation of the number of neighbors. 




0.1 0.2 0.3 0.4 0.5 

AR/R 

Fig. 10. Area under the first peak of the pair correlations of 
the three contact types. Symbols follow fig. O 



One can directly infer that up to intermediate dispersi- 
ties, where the critical packing fraction reaches its lowest 
value, small particles increase their probability of being 
around other particles, of both types (top and bottom pan- 
els) with a more pronounced effect at / + = 0.800, while 
the probability for large-large particle contacts reaches its 



lowest value (middle panel). In addition to that, the large- 
large contacts barely change with number fraction and its 
minimum value occurs at a dispersity value close to those 
of the minimum jamming density, fig. [3l 

These results, along with the RDFs in figs. [6] and [71 im- 
ply that, less efficient packing correspond to more small 
particle clusters and less large particle ones. A possible 
explanation can be given by the fact that the initial com- 
pression is provided solely by large particles and the pack- 
ing relaxation should take the available space provided 
by small-small neighbors, i.e., large particles should push 
small ones in order to decrease the potential energy. This 
will deform the initial structure and small-small contacts 
will be formed in a distinct structure that the initial one. 
Since this structure is the most dense possible, the critical 
packing will occur at a lower density. Also, the larger the 
number of small particle contacts, the larger is the space 
available for structure rearrangements. Therefore, less de- 
compression cycles will be needed to reach the minimum 
energy. 

For polydisperse packings, the jammed structure shows 
a completely distinct scenario. Since there are several par- 
ticle sizes, the chance for the formation of crystalline re- 
gions during the packing relaxation is very low. There- 
fore, the packing structure should be strongly amorphous, 
as shown in fig. [11] at large AR/Rq. Moreover, this amor- 
phization seems to be a continuous process, since the peaks 
at low dispersity become smoother for larger AR/Rq until 
they merge and, eventually, disappear. 
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Fig. 11. Pair correlation of multidisperse packings for disper- 
sity values (from top to bottom): AR/Ro = 0.500, 0.200, 0.120, 
0.050 and 0.010. Curves are shifted for clarity. 

Since the results for the RDFs implied that the small 
particle contacts introduced disorder, one can find a corre- 
lation between the CPF and the contacts orientational or- 
der. Fig.[T2]has the average value, over particles and runs, 
measured at the end of the decompression, of the orien- 
tational order parameter related to the triangular lattice, 
eq. (|3]). It should be noticed that, given the definition of 
first neighbors followed here, this order parameter does 
not contain any contribution from rattlers. 

First of all, one see that, at low dispersity, this order 
parameter decreases continuously. This is consistent with 
the arguments given earlier, that small particle contacts 
introduce disorder in the system and it packs less effi- 
ciently. Moreover, the fast increase of the critical packing 
fraction at large dispersities and number fractions can be 
seen to correlate with a fast increase in the orientational 
order, also consistent with the appearance of large crys- 
talline regions. On the other hand, the smooth increase 
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Fig. 12. Orientational order parameter for all number frac- 
tions and dispersity values. Symbols and colors are as in fig. 
rj /+ = 0.400 (circles), /+ = 0.500 (squares), /+ = 0.600 (di- 
amonds), /+ = 0.700 (triangles), /+ = 0.800 (left triangles), 
polydisperse (inverted triangles). 

in (</>) with dispersity is not accompanied with a corre- 
sponding increase in the order parameter. This implies 
that the inversion in the CPF curves are due to pockets 
of large particle clusters, as seen in fig. [TU] (middle panel), 
regardless of the overall decrease in order. Large parti- 
cles clusters have a high value of \(&)\, and, given the 
smooth increase in Ni,l(t) with dispersity, one may infer 
that their number is small, since the whole system has a 
low value of the order parameter. Only when large particle 
clusters are formed, the global order increases, providing 
a more dense packing. 

In the polydisperse case, the orientational order pa- 
rameter continuously decreases with increasing dispersity. 
Its value is higher than for the bidisperse cases up to 
AR/Ro = 0.200. This is surprising for one would ex- 
pect that more particle sizes would lead to less order 
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(from eqs. ([5] and (O) one can see that the size disper- 
sion of a polydisperse packing is proportion to AR) . Given 
that there is no monodisperse limit for this packing when 
AR/ Rq — > 1, one can imagine that the order parameter for 
the f + = 0.400 and 0.500 cases will eventually be higher 
than the polydisperse one at some dispersity beyond 0.500. 

4 Summary and conclusions 

It was presented a numerical study on the jamming prop- 
erties and structure of a two dimensional packing of elas- 
tic disks, for bi- and polydisperse cases. The attention was 
focused on the value of the maximal packing fraction for 
which the compression energy is zero. This was measured 
through numerical decompression experiments of a disor- 
dered packing, initially arranged in a crystal (triangular 
lattice) structure. The critical packing fraction (CPF) was 
measured at the end of the decompression and was shown 
as a function of the dispersity degree, AR/Rq. Also, for 
bidisperse packings, the CPF was also studied as a func- 
tion of the number fraction of large disks. 

The general trend of the CPF is initially decreasing 
up to a minimum value, and then increasing with disper- 
sity, fig. [3J The lowest CPF value is close to the RCP 
value, but obtained for a number fraction of /+ = 0.400, 
instead of the original value obtained at / + — 0.500 [2]. 
The distinct CPF behavior with AR/R observed in [TUl 
[T7] is due to the fact that, in the present case, one starts 
from the most possible dense packing and, therefore, the 
introduction of disorder, through the quench and inter- 
nal structure rearrangement, will certainly lead to a lower 



a novel route to unjamming in soft disk packings 13 

jamming packing fraction, since a more efficient packing 
can be achieved with a increase in local order |24) . 

At low dispersity, the system behaves approximately 
as the monodisperse crystal, as seen in the results for the 
average number of decompression steps, fig. HI The depar- 
ture from the monodisperse regime can be attributed to 
an increase in the number of small-small particle contacts, 
ng.EOl 

The structure reveals that, for low dispersity, the de- 
compressed packing has significant order, as revealed by 
the long range behavior of gss( r ), 9ll(j~) and gsL{r). The 
packing structure is mostly disordered at intermediate dis- 
persities, and at the highest dispersity, large-large particle 
contacts bear most of the translational order in the sys- 
tem, which forms a crystal with defects, with the few small 
particles scattered between the scarce space available be- 
tween the large particle contacts. 

For polydisperse packings the long range order is com- 
pletely absent for large dispersities, since a disk size is 
chosen from a uniform distribution of values in the range 
[—AR, AR], for larger AR, there will a very broad distri- 
bution of particle sizes. 

The data for the local orientational order gives a sim- 
ilar picture. For larger number fractions, the packings are 
more ordered locally compared to lower number fraction 
cases. Also, for low dispersities, the packing is more or- 
dered, as expected since it behaves like the monodisperse 
one. Above the low dispersity range, the local orientational 
order decreases with AR/ Rq, regardless of the smooth in- 
crease in the CPF. However, the fast increase of the CPF 
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at large number fraction and dispersities is strongly corre- 
lated with a fast increase in this order parameter, implying 
that such packings are close to a crystal. 
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